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■ Abstract 

o 

Model of a passive scalar field advected by the compressible Gaussian strongly aniso- 
tropic velocity field with the covariance oc 8(t — t')l x — x'| 2e is studied by using the field 



theoretic renormalization group and the operator product expansion. The inertial-range 
stability of the corresponding scaling regime is established. The anomalous scaling of the 
single-time structure functions is studied and the corresponding anomalous exponents are 



calculated. Their dependence on the compressibility parameter and anisotropy parameters 
is analyzed. It is shown that, as in the isotropic case, the presence of compressibility leads 
to the decreasing of the critical dimensions of the important composite operators, i.e., the 
anomalous scaling is more pronounced in the compressible systems. All calculations are 
done to the first order in e. 
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1 Introduction 



During the last two decades the so-called "rapid change model" of a passively advected scalar by 
a self-similar Gaussian 5— correlated in time velocity field introduced by Kraichnan [1] and num- 
ber of its extensions have played the central role in the theoretical investigation of intermittency 
and anomalous scaling, the problems which stay in the center of attention in the framework of 
the inertial range investigation of fully developed turbulence [21 [3], i.e., in the range character- 
ized by the scales which are far away from the largest scales at which energy is pumping into the 
system and, at the same time, far away form the smallest scales which are related to the dissi- 
pation processes. There are at least three reasons for this interest: First of all, it is well-known 
from experimental and also theoretical studies that the deviations from the statements of the 
famous classical Kolmogorov-Obukhov phenomenological theory (see, e.g., [21 HIE]) is surpris- 
ingly more noticeable and visible for the simpler models of passively advected scalar quantity 
(scalar field) than for the velocity field itself (see, e.g., [SJ [TJ, [HI El HDl EH]); second, at the same 
time, the problem of passive advection of scalar field (as well as vector field) is much easier 
from theoretical point of view than the original problem of anomalous scaling in the framework 
of Navier-Stokes velocity field, and, in the end, third, even very simplified models with given 
Gaussian statistics of velocity field lead to the anomalous behavior which describe many fea- 
tures of the real turbulent advection, see, e.g., Refs. [II El El [IQl [III [121 [131 [III [ISl [111 [13 [IE] and 
references cited therein. As was already mentioned, the crucial role in these studies was played 
(and is still played) by the aforementioned rapid change model of a passive scalar advection, 
where, for the first time, the systematic analysis of the corresponding anomalous exponents was 
done on the microscopic level. For example, within the so-called "zero-mode approach" to the 
rapid change model [12] (see also survey paper [19]) the anomalous exponents are found from 
the homogenous solutions (zero modes) of the closed equations for the single-time correlations. 

On the other hand, one of the most effective approach for studying self-similar scaling 
behavior is the method of the field theoretic renormalization group (RG) [20, 2lJ. It can be 
also used in the theory of fully developed turbulence [221 [23] an d related problems, e.g., in the 
problem of a passive scalar advection by the turbulent environment [23] (see also Refs. [21l l25l 
ES] for details). 

In Refs. [27] the field theoretic RG and operator-product expansion (OPE) were used in 
the systematic investigation of the rapid-change model of a passive scalar. It was shown that 
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within the field theoretic approach the anomalous scaling is related to the very existence of 
so-called "dangerous" composite operators with negative critical dimensions in OPE (see, e.g., 
Refs. [21| [26] for details). In the subsequent papers [28] the anomalous exponents of the model 
were calculated within the e expansion up to order e 3 (three-loop approximation). Here e is 
a parameter which describes a given equal-time pair correlation function of the velocity field 
(see next section). Important advantages of the RG approach are its universality and calcu- 
lational efficiency: a regular systematic perturbation expansion for the anomalous exponents 
was constructed, similar to the well-known e-expansion in the theory of phase transitions. 

Besides various generalizations of the Kraichnan model towards more realistic ones, namely, 
models with inclusion of small scale anisotropy [29], compressibility [301 13T]. an d finite correla- 
tion time of the velocity field [321 [331 El] were studied by the field theoretic approach. General 
conclusion of all these investigations is that the anomalous scaling, which is the most intriguing 
and important feature of the Kraichnan rapid change model, remains valid for all generalized 
models. 

In Ref. [29] the field theoretic RG and OPE were applied to the rapid change model of passive 
scalar advected by Gaussian strongly anisotropic velocity field where the anomalous exponents 
of the structure functions were calculated to the first order in e expansion. It was shown that 
in the presence of small-scale anisotropy the corresponding exponents are nonuniversal, i.e., 
they are functions of the anisotropy parameters, and they form the hierarchy with the leading 
exponent related to the most "isotropic" operator. The importance of these investigations is 
related to the question of the influence of anisotropy on inertial-range behavior of passively 
advected fields as well as the velocity field itself [ElQSl[32l[33l[3g (see also 

the survey paper [H] and references cited therein, as well as recent astrophysical investigations, 
e.g, in Refs. [121 S3]). On one hand, it was shown that for the even structure (or correlation) 
functions the exponents which describe the inertial-range scaling exhibit universality and they 
are ordered hierarchically in respect to degree of anisotropy with leading contribution given by 
the exponent from the isotropic shell but, on the other hand, the survival of the anisotropy in 
the inertial-range is demonstrated by the behavior of the odd structure functions, namely, the 
so-called skewness factor decreases down the scales slower than expected earlier in accordance 
with the classical Kolmogorov-Obukhov theory. 

On the other hand, in Ref. [31] the influence of compressibility and large-scale anisotropy on 
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the anomalous scaling behavior was studied in the aforementioned model in the order e 2 in the 
e expansion, and anomalous exponents of higher-order correlation functions were calculated 
as functions of the parameter of compressibility. It was shown that the exponents exhibit 
hierarchy related to the degree of anisotropy. Again, the existence of small-scale anisotropy 
effects were demonstrated by the odd dimensionless ratios of correlation functions: the skewness 
and hyperskewness factors, and it was shown that the persistent of small-scale anisotropy is 
more pronounced for larger values of the compressibility parameter a. From this point of view, 
compressible systems are very interesting to be studied. 

In present paper, we shall continue in the investigation of the model, namely, for the first 
time, we shall study the model with compressible velocity field together with assumption about 
its small-scale uniaxial anisotropy in the first order in e, i.e., we shall extend the model studied 
in Ref. [29J to the compressible case. In this situation, in general, two types of diffusion- 
advection problems of scalar particles exist in the compressible case which are identical in the 
incompressible case. In what follows, we shall study only one of them, namely, the advection 
of so-called "tracer" (see, e.g., Ref. [31]), i.e., for example, concentration of a scalar impurity, 
temperature, entropy, etc. Details will be shown in the next section. 

First of all we shall establish stability of the scaling regime of the model and coordinates 
of the corresponding infrared (IR) stable fixed point will be found analytically as functions of 
the compressibility and anisotropy parameters. These results will be then used in the analysis 
of the asymptotic behavior of the single-time structure functions of a passively advected scalar 
field. 

The paper is organized as follows: In Sec. [2] we present definition of the model and introduce 
the compressibility and small-scale uniaxial anisotropy to the given pair correlation function 
of the velocity field. In Sec. [3] we give the field theoretic formulation of the original stochastic 
problem and discuss the corresponding diagrammatic technique. The analysis of the ultraviolet 
(UV) divergences of the model is given, the multiplicative renormalizability of the model is 
established, and the renormalization constants are calculated in one-loop approximation. In 
Sec. H] we analyze infrared (IR) asymptotic behavior of the model which is governed by the IR 
stable fixed point. Explicit expressions for the coordinates of the IR fixed point are found. In 
Sec. the renormalization of needed composite operators is done and their critical dimensions 
are found as functions of parameters of the model. In Sec. [6] discussion of results is present. 
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2 The model of advection of passive "tracer" 

We shall consider the problem of the advection of a passive "tracer" 9{x) = 6{t,x) which is 
described by the following stochastic equation 

d t e = v A9-(v i d i )d + f, (1) 

where dt = d/dt, di = d/dxi, A = d 2 is the Laplace operator, vq is the coefficient of molecular 
diffusivity (hereafter all parameters with a subscript denote bare parameters of unrenormal- 
ized theory; see below), Vi = Vi(x) is the i-th component of the compressible velocity field v(x), 
and / = f(x) is a Gaussian random noise with zero mean and correlation function 

D f = (f(x)f(x')} = 5(t - t')C(r/L), r = x - x', (2) 

where parentheses (...) hereafter denote average over corresponding statistical ensemble. The 
noise defined in Eq. (T5]) maintains the steady-state of the system but the concrete form of the 
correlator will not be essential in what follows. The only condition which must be fulfilled by 
the function C(r/L) is that it must decrease rapidly for r = |r| ^> L, where L denotes an 
integral scale related to the stirring. 

As was already mentioned in Introduction, another type of the diffusion-advection problem 
exists in the case when compressibility of the velocity field is supposed which is given by the 
following more general stochastic equation [311 IH] 

d t 9 = u A9 - diiyiO) + f. (3) 

It describes the passive advection of a density field but, in what follows, we shall not study 
this problem. In the case of incompressible velocity field (it is given mathematically by the 
divergence-free condition <9jt>j = 0) both models are equivalent. 

In real problems it is traditionally assumed that the velocity field v(x) satisfies stochastic 
Navier-Stokes equation [21]. In spite of, in what follows, we shall suppose that the velocity field 
obeys a Gaussian distribution with zero mean and two-point correlator 

(v t (x) Vj (x')) = Dl{x ]X >) = D 5(t-f) J |^-- exp[ik(x-xO], (4) 

i.e, we shall work in the framework of the so-called rapid-change model [H [T21 [E3 [HI [23 I2H 
[29| [30l HH HHl H7] . Here, d denotes the dimension of the x space, and D is an amplitude factor 
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related to the coupling constant g of the model (expansion parameter in the perturbation 
theory, see next Section) by the relation D /u = g ~ A 2e , where A is the characteristic UV 
momentum scale. The parameter of the energy spectrum of the velocity field < e < 1 is taken 
in such a way that its "Kolmogorov" value (the value which corresponds to the Kolmogorov 
scaling of the velocity correlation function in developed turbulence) is e = 2/3, and 1/m is 
another integral scale. In general, the scale 1/m may be different from the integral scale L 
introduced in Eq. fl2]) but, in accordance with Ref. [29J, we suppose that 1/m ~ L. 

In the incompressible isotropic case the second-rank tensor Rij(k) in Eq. (J4j) has the simple 
form of the ordinary transverse projector -Rjj(k) = i^j(k) = 5^ — kikj/k 2 , where k = |k|. 
This tensor is changed when one incorporates small-scale anisotropy or compressibility. Let us 
briefly discuss these questions. 

First of all, in the case of incompressible anisotropic case a new second-rank tensor must 
be again a transverse operator because the velocity field is still divergence-free (c^fj = 0). The 
simplest way how to introduce small-scale uniaxial anisotropy is to take the operator Rij(k) in 
the following way [29] 

%(k) = ^1 + ai ^j^j i^(k) + a 2 n a n,P ia (k)P J -,(k) , (5) 

where -Py(k) is the usual transverse projection operator (as defined above), the unit vector n 
determines the distinguished direction, and a\, a2 are parameters characterizing the anisotropy. 
The positive definiteness of the correlation function (j3J) imposes the following restrictions on 
their values: ot\ , ol% > —1. The operator (111) is a special case of the general transverse structure 
that possesses uniaxial anisotropy: 

i2y(k) = a(V0^i(k) + 6(V')n a n l P fa (k)P i ,(k) , (6) 

where ip denotes the angle between the vectors n and k (n ■ k = kcosip). Using Gegenbauer 
polynomials [H] the scalar functions in representation ([6]) may be expressed in the form 

oo oo 

a W = aiF>2i ( COS ^) ' K4>) =Y1 b i P 2i (COS *0) . 
i=0 i=0 

It was shown in Ref. [29J that all main features of the general model with the anisotropy 
structure represented by Eq. are included in the simplified model with the special form of 
the transverse operator given by Eq. (151) . 
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Second, in the case of the compressible isotropic velocity field the second-rank tensor Rij(k) 
in Eq. (j3J) is not longer a transverse projector but it is rather a combination of the transverse 
projector -Pij(k) = 8y — kikj/k 2 and of the longitudinal projector Qij(k) = kikj/k 2 as a result 
of the fact that the velocity field is not longer solenoidal one. Thus, in this case, the correlator 
(Tj0) contains the tensor structure of the following general form 

R ij (k) = P ij (k)+aQ ij (k) (7) 

where a > is a free parameter of the compressibility. The value a = corresponds to 
the divergence-free (incompressible) advecting velocity field. It was used in Refs. [30], [31] for 
determination of the influence of the compressibility on the anomalous scaling of the correlation 
functions of scalar density, as well as tracer fields in two-loop level. 

In what follows, we shall study the combined effects given by the small-scale anisotropy 
and compressibility and our aim will be to study possible deviations from the conclusions given 
in Ref. [29], where the rapid-change model of passively advected scalar field with small scale 
uniaxial anisotropy was studied. To do this, it is necessary to introduce into the velocity 
field correlator (J4]) corresponding second-rank tensor which will have needed properties. The 
simplest way how to do this is to add longitudinal projector Qij(k) = kikj/k 2 to the uniaxial 
anisotropic transverse tensor structure given in Eq. Thus, the result tensor structure which 
will be used in our investigations is defined as 

fn ■ k) 2 

i^-(k) = P y (k) + aQtf(k) + V / P i3 -(k) + a 2 n 8 n,P i8 (k)P j7 (k) . (8) 

It means that we have introduced the compressible term to the isotropic component of the 
axially anisotropic tensor structure (0) only. Therefore, the contributions of compressibility 
and anisotropy are given by a simple sum of the corresponding terms. Nevertheless, as we shall 
see, the combined effects of anisotropy and compressibility on the results will not be a simple 
sum of them. 



3 Field Theoretic Formulation of the Model, UV Renor- 
malization and RG analysis 

According to the well-known general theorem (see, e.g., Refs. [20] |2T] ) the stochastic problem 
(JTJ) and ([2]) is equivalent to the field-theoretic model of the set of three fields $ = {8, 6', v} with 

7 



the following action functional 

= ~\J dt 1 d d x 1 dt 2 d d x 2 v i (t 1 ,x 1 )lD^t 1 ,x 1 -,t 2 ,x 2 )}- 1 v j ( y t 2 ,x 2 ) 
+^ J dt x rf d Xl dt 2 d d x 2 9'{t u *i)D f (t u x l5 1 2 , x 2 )0'(t 2 , x 2 ) 
+ /" dt d d x 0' [-<9 t - u^j + i/ A] 9, (9) 

where 0' is an auxiliary field, and all summations over the vector indices are implied. The 
second and the third integral in Eq. ([9]) correspond to the Martin-Siggia-Rose action [49J for 
the stochastic problem (JTJ), (j2J) at fixed velocity field v, and the first integral describes the 
Gaussian averaging over v defined by the correlator D" in Eq. (j3J) with tensor Rij(k) given in 
Eq.©. 

Action ([9]) is given in a form convenient for application of the field theoretic perturbation 
analysis with the standard Feynman diagrammatic technique. From the quadratic part of the 
action one obtains the matrix of bare propagators. The wave-number-frequency representation 
of propagators of the fields 6 and 9' is 

C(k) (ii) 



(9'6') = (12) 

where C(k) is the Fourier transform of the function C(r/L) form Eq. On the other hand, the 
bare propagator of the velocity field (vv)q is defined by Eq. (J3J) with the tensor structure given 
by Eq. (jSJ). In what follows, we shall need only the propagators (99') o and (vv)o- Their graphical 
representation is shown in Fig.[TJ The interaction in the model is given by the nonlinear term 
—9'(vidi)9 = 9 r VjVj9 with the vertex factor which in the wave-number-frequency representation 
has the form (the momentum flows into the vertex via the scalar field 9): 

Vj = -iky (13) 

Its graphical representation is given in Fig.[TJ 

Standard power counting [20, |21] leads to the identification of correlation functions with 
superficial UV divergences. In the framework of the rapid-change passive advection models 
detail analysis of this question was done, e.g., in Ref. [29], where it was shown that the only 
one-particle-irreducible (1PI) Green function which possesses superficial UV-divergences is the 
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Figure 1: (Left) The graphical representation of needed propagators of the model. (Right) The 
interaction vertex of the model (wave-number-frequency representation). The momentum flows 
into the vertex via the scalar field 9. 



function Tgig = (9'6)i^i r (within the rapid-change model the situation is unchanged when com- 
pressibility of the system is assumed) and, in the isotropic case, this Green function leads only 
to the renormalization of the term i>q9'A9 of action (Q and the corresponding UV divergences 
may be fully absorbed in the adequate redefinition of the existing parameters go, u . Thus, all 
correlation functions calculated in terms of the renormalized parameters g, v are UV finite. 

The situation becomes, however, more complicated when anisotropy is introduced. It is 
related to the fact that in this case the 1PI Green function Tg/g contains divergences corre- 
sponding to the structure ^'(n • d) 2 9 (the only possible anisotropic structure ~ A = d 2 ) which 
is not present in original unrenormalized action (Q. It leads to the non-renormalizability of the 
model in the anisotropic case. To make the model multiplicatively renormalizable it is necessary 
to extend original action by including needed term with corresponding new parameter. As 
a result the extended model is described by the following action 



where a new unrenormalized parameter xo has been introduced. As was pointed out in Ref. [29J 
the stability of the system requires the positivity of the total viscous contribution z/ /c 2 + Xo I/ o( n " 
k) 2 , i.e., the inequality Xo > — 1 must be fulfilled. This modification leads, of course, also to 
the modification of the corresponding isotropic propagators of the fields 9, 9' given in Eqs. fllO|) 
and (ITT]) which are now defined as (see also Ref. [29J) 




(14) 



After this modifications the model defined by action ([HI) becomes multiplicatively renor- 
malizable and the standard RG analysis can be now applied. The corresponding renormalized 
action has the form 

S R {$) = --J dt 1 c? d x 1 dt 2 d d x 2 ^(t 1 ,x 1 )[ J D^.(t 1 ,x 1 ;t 2 ,x 2 )]- 1 ^(t 2 ,x 2 ) 

+ i J dh d d Xl dt 2 d d x 2 XijD^Cti, xi; t 2 , x 2 )0'(t 2j x 2 ) (17) 
+ J dtd d *9' [-d t - vA + uZ x A + yyZ 2 ip. ■ d) 2 ] 0, 

where Zi and Z 2 are the renormalization constants (they absorb the UV divergent parts of the 
1PI function Tqiq). It is equivalent to the multiplicative renormalization of the bare parameters 
0o, i/ , and Xo, namely 

i/ = 5-o = fi7i 2 %, xo = XZ X , (18) 

where 0, u, and x are renormalized counterparts of the corresponding bare parameters, and /i is 
a scale setting parameter or the reference mass (it has the same canonical dimension as the wave 
number). In what follows, we shall work in minimal subtraction scheme (MS), therefore, in one- 
loop approximation, the renormalization constants Z have the form 1 + A(g, a, x, «i, ot 2 , d)/e. 
Thus, A is a function of dimensionless parameters but it is independent of e. 

By comparison of the corresponding terms in action (IT7|) with definitions of the renormal- 
ization constants Z for parameters ([TBI , one obtains the following relations between renormal- 
ization constants: 

Z v = Zi, Z x = Z 2 Z 1 , Z g = Z x , (19) 

where last relation in Eq. (I19p is a consequence of the fact that D defined in Eq. (J3J) is not 
renormalized, i.e., D = g§v§ = gv\£ e . 

Standardly, the formulation through the action functional (Q (or through (fTTl) in the 
anisotropic case) replaces the statistical averages of random quantities in the stochastic problem 
defined by Eqs. flTJ and (|2J) with equivalent functional averages with weight expS'($). Gener- 
ating functionals of total Green functions G(A) and connected Green functions W(A) are then 
defined by the functional integral 

G(A) = e w{A ^ = Jv$ e s W +A * } (20) 

where A(x) = {A e , A e ' , A v } represents a set of arbitrary sources for the set of fields $, £>$ = 



10 



T>9T>9'T>v denotes the measure of functional integration, and linear form A$ is defined as 

A$ = J dx[A\x)6(x) + A 6 '(x)B'{x) + A v t {x)vi{x)}. (21) 

Let us continue with renormalization of the model. The relation S(9, 9', v, eo) = S R (8, 6', v, e, //), 
where eo stands for the complete set of bare parameters and e stands for renormalized one, leads 
to the relation W(A, eo) = W R (A, e, fi) for the generating functional of connected Green func- 
tions. By application of the operator = \xdp_ at fixed eo on both sides of the latest equation 
one obtains the basic RG differential equation 

V RG W R (A,e,n)=0, (22) 

where T)rg represents operation T>^ written in the renormalized variables. Its explicit form is 

V RG = V fl + p g (g, x)d g + P x {g, x)d u - l»(g, x)V», (23) 

where we denote D x = xd x for any variable x and the RG functions (the (3 and 7 functions) 
are given by well-known definitions (see, e.g., Refs. [201 [21]) an d, in our case, by using relations 
ffT9]) for renormalization constants, they have the following form 

7< = V^lnZi, i = l,2 (24) 
f3 g = f> fl g = g(-2e + ll ), (25) 
Px = ^X = x(7i-7 2 ). (26) 

The renormalization constants are determined by the requirement that the one-particle 
irreducible Green function (#'0)i_j r must be UV finite when is written in renormalized variables. 
In our case it means that they have no singularities in the limit e — > 0. The one-particle 
irreducible Green function (9'9)i^ ir is related to the self-energy operator by the Dyson 
equation 

(d'd)i-ir = -iu + vqV 2 + ^oXo(n • p) 2 - T, e > e {u,p). (27) 

Thus Z\ t 2 are found from the requirement that the UV divergences are canceled in Eq. (I27p 
after substitutions z/ = vZ\ and z/oXo = vx^i- This determines Z\^ up to an UV finite 
contributions, which are fixed by the choice of the renormalization scheme. In the MS scheme, 
as was mentioned already above, all the renormalization constants have the form: 1 + pole in 
e. In one-loop approximation the self-energy operator Eg/g is represented by the corresponding 
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£(9'<9 — 



Figure 2: The only diagram which contribute to the self-energy operator Yle'o- 

one-particle irreducible diagram which is shown in Fig. [2] but it must be stressed that, at the 
same time, it is an exact result because in rapid-change model all higher-loop diagrams contain 
at least one closed loop which is built on by retarded or advanced propagators only, thus they 
are automatically equal to zero. 

The diagram in Fig.[2]has the following analytical representation 

^=TsMK£n<^ + <■•■»>'*>■ < 28) 

where Sd = 2ir d ^ 2 /T(d/2) denotes the d-dimensional sphere, and the functions A\ and A2 have 
the following simple explicit form 

Al ~ l(dT2) ' (29) 

_ -2(ai + a 2 ) + d 2 a 2 , , 

A " ~ dJdTv ' (30) 

In the end, the renormalization constants Z\ and Z 2 are given as follows 

(3D 

Z 2 = l- 9 - 7 %^. (32) 

Because of the second rank tensor (jSJ) is taken in the simple form of the sum of the uniaxial 
anisotropic transverse projector and the longitudinal projector then, as a result of this fact, the 
self-energy operator T, g / d is also a simple sum of the results obtained in Refs. [29] and [31]. Using 
the definition of anomalous dimensions 71^ in Eq. (|24"I) one comes to the following expressions 

7i = 77^1 , 72 = 77-^2, (33) 
2 2x 

where we denote g = gSd/(2n) d , and A 1<2 are given in Eqs. (1291 and (l3"Uj) . 

Let us note once more that the expressions fl311 - fl33l are exact as a result of non-existence 
of the higher-loop corrections. 
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4 Fixed point and scaling regime 

Possible scaling regimes of a renormalizable model are directly given by the IR stable fixed 
points of the corresponding system of RG equations [201 121] • The fixed point of the RG 
equations is defined by /3-functions, namely, by requirement of their vanishing. In our model, 
the coordinates g*, x* of a fixed point are found from the system of two equations 

Pg(.9*, X*) = i3 x (g*, X*) = 0. (34) 

The beta functions (3 g and (3 X are defined in Eqs. (j25|) . and (1261) . To investigate the IR stability 
of a fixed point it is enough to analyze the eigenvalues of the matrix Q of the first derivatives: 



(35) 



d/3 g /dg d(3 g /d X 

The IR asymptotic behavior is governed by the IR stable fixed points, i.e., those for which 
both eigenvalues are positive. Using the explicit expressions given in Eq. (!33l) together with 
the definitions of f3 functions in Eqs. ( 1251) and ( 1261) leads to the explicit expressions for the 
coordinates of the non-trivial IR stable fixed point 

4d(d + 2)e 

9 * "~ {d + 2)(d-l + a) + a l {d+l) + a 2 ) 1 } 

-2(ai + a 2 ) + d 2 a 2 

X * ~ (d + 2)(d-l + a) + a 1 (d+l) + a 2 ' [ } 

In the incompressible limit a — > one comes to the results of Ref . [2H] , namely 



-a = Hd + 2)e 

y * ' {d + 2){d-l) + a x (d+l) + ot2 1 J 

-2(ai + a 2 ) + rf 2 a 2 
X * (d + 2)(d-l) + ai (d+l) + a 2 ' 

On the other hand, in isotropic limit «i i2 — >■ one has (see, e.g., Ref. [33] ) 

c = 4d(d + 2) £ 
y * (d + 2)(d-l + a) 1 J 

Thus, as one can see, our definition of the second rank tensor ([8]) leads to the simple relation 
among g*,J%, and g%, namely 

111 

— = — + -. 41 

g* m m 
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The values of anomalous dimensions 71 and 72 are found exactly at fixed point and are 
defined as follows 

it = 7 2 * = 2e, (42) 

where 7^ = 71,2(5-*, X*)- 

The matrix of the first derivatives taken at the fixed point has simple diagonal form 

( 2e \ 

= , (43) 

\0 2s) 

i.e., the fixed point is IR stable if e > 0. This is also the only condition to have g* > (together, 
of course, with earlier discussed physical assumptions, namely: a.\ , 2 > — 1 and a > 0, see Sec.|2]). 
The physical condition > — 1 is also fulfilled without further restrictions on the parameter 
space. 

The issue of interest are especially multiplicatively renormalizable equal-time two-point 
quantities G(r) (see, e.g., [2D]). The example of such quantity are the equal-time structure 
functions 

S N (r) = ([9(t,x)-e(t,x')} N ), r=|x-x'| (44) 

in the inertial range specified by the inequalities / ~l/A<<r<<L = l/m(Zisan internal 
length). Here parentheses (■ • •) mean functional average over fields <£> = {8, 8', v} with weight 
exp S R ($). 

First let us describe briefly IR scaling behavior in general on the example of an equal-time 
function G(r) which is multiplicatively renormalizable. The IR scaling behavior of the function 
G(r) (for r/l ^> 1 and any fixed r/L) 

G{r) ~ uf G r da {r/l)- AG R{r/L) (45) 



is related to the existence of IR stable fixed points of the RG equations (see above). In Eq. fj45|) 
cIq and dc are corresponding canonical dimensions of the function G: ($q is the frequency 
dimension, and do is the total canonical dimension. They are related by the relation do = 
d@ + 2d,Q, where dg, is corresponding momentum canonical dimension. The existence of two 
independent canonical dimensions is related to the fact that our model belongs among two-scale 
dynamical models (details see, e.g., in Refs. [HI [26]). In Eq. (|45l) R{r / L) is so-called scaling 
function which cannot be determined by RG equation (see, e.g., [21]), and A G is the critical 
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dimension denned as 

A G = d k G + A^ G + 7* . (46) 

Here 7q is the fixed point value of the anomalous dimension 7^ = fid^lnZc, where Zq is 
renormalization constant of multiplicatively renormalizable quantity G, i.e., G = ZqG r |33j, 
and A w = 2 — 7* = 2 — 7J 1 is the critical dimension of frequency with 7^ = 2e as it is shown in 
Eq.(H2]). 

Now, let us apply the above discussion to the inertial-range analysis of the equal-time 
structure functions as defined in Eq. f|44l . It is well-known that, in the isotropic case, the odd 
functions S2 n +i vanish, while for S2 n simple dimensional considerations give 

S 2 n{r) = v n r 2n R 2n (r/l,r/L), (47) 

where Rm are some functions of dimensionless variables. First, the multiplicative renormal- 
izability of the model leads to the existence of differential RG equations for these structure 
functions and their asymptotic behavior for r/l » 1 and any fixed r/L is given by IR stable 
fixed point of the RG equations and the structure functions can be written in the following 
form 

5 2 „(r) = v n r 2n (r/l)-^R 2n (r/L), r/l » 1 (48) 

with unknown scaling functions R2 n {r/L). In the theory of critical phenomena [201 121] the 
quantity 7„ is known as "anomalous dimension" and A[S2„] = — 2n + 7„ is termed the "critical 
dimension" (see above), where — 2n is the corresponding "canonical dimension". In our case, 
7„ = 2ne [21], i.e., representation (j4~8l) implies the existence of a scaling in the IR region 
(r/l » 1, r/L fixed) with definite critical dimensions of all "IR relevant" parameters, AfSy = 
— 2n(l — e), A r = —1, A L -i = 1 and fixed "irrelevant" parameters uq and / regardless of the 
form of the functions R2 n {r/L). 

The second stage of RG analysis is associated with the investigation of small r/L behavior 
of the functions R2 n (r/L) in Eq. (148]) using the OPE. It shows that, in the limit r/L — >• 0, the 
functions R2 n (r/L) have the asymptotic form 

R 2n (r/L) = Y / C F (r/L)(r/L) A ", (49) 

F 

where Cp are some coefficients regular in r/L and the summation is implied over certain 
renormalized composite operators F with critical dimensions A n [21J. In the case under consid- 
eration, the leading operators F have the form F n = (di6di6) n . In Sec.[5]we shall consider them 
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Figure 3: Graphical representation of the one-loop correction to the composite operator T N 
(details see in Ref. [29]). 

in detail where the complete calculation of the critical dimensions of the composite operators 
F n will be present for arbitrary values of n, d, a and a 12 - 

5 Operator product expansion, Critical dimensions of 
composite operators, and Anomalous scaling 

5.1 Operator product expansion 

Let us now study the behavior of the scaling function in Eq. fj45|) . According to the OPE 
J20 , EH ESI EE] , the equal-time product Fi(x')F 2 (x") of two renormalized composite operators 
□ at x = (x' + x")/2 = const and r = x' — x" — > can be written in the following form 

F 1 (x')F 2 (x") =J2C Fi (r)Fi(x,t), (50) 

i 

where summation is taken over all possible renormalized local composite operators Fi allowed 

by symmetry with definite critical dimensions A p. , and the functions Cp i are the corresponding 

Wilson coefficients regular in L~ 2 . The renormalized correlation function (F\{x')F 2 (x")) can be 

now found by averaging Eq. fl50|) with the weight exp S R with S R from Eq. (1171) . The quantities 

(Fi) appear on the right-hand side, and their asymptotic behavior in the limit L^ 1 — ^ is then 

found from the corresponding RG equations and has the form (Fi) oc L~ Ap i . 

From the OPE (!50l) one can find that the scaling function R(r/L) in the representation (|45|) 

2 By definition we use the term "composite operator" for any local monomial or polynomial constructed from 
primary fields and their derivatives at a single point x = (t,x). Constructions 9 n (x) and [di9(x)did(x)] n are 
typical examples 
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Figure 4: Dependence of the critical dimension A [2, 2]/e on anisotropy parameters a,\ (a 2 = 0) 
and ot 2 {at\ = 0) for different values of the compressibility parameter a. 

for the correlation function F\{x')F 2 (x") has the form 

R(r/L)=J2C Fl (r/L)^, (51) 

i 

where the coefficients C Fi are regular in (r/L) 2 . 

The principal feature of the turbulence models is the existence of the so-called " dangerous" 
operators with negative critical dimensions [2T1 1251 1261 [271 [30] . Their presence in the OPE 
determines the IR behavior of the scaling functions and leads to their singular dependence on 
L when r/L —>■ 0. Therefore, the turbulence models are crucially different from the models of 
critical phenomena, where the leading contribution to the representation f[5"Tj) is given by the 
simplest operator F — 1 with the dimension Ap = 0, and the other operators determine only 
the corrections that vanish for r/L — > 0. 

If the spectrum of the dimensions for a given scaling function is bounded from below, 
the leading term of its behavior for r/L —>■ is given by the minimal dimension. It will be 
our case and for the investigation of the anomalous scaling of the structure functions (|44|) 
the leading contribution of the Taylor expansion is given by the tensor composite operators 
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Figure 5: (Left) Dependence of the critical dimension A[2, 2]/e on anisotropy parameters ot\ = 
cli for different values of the compressibility parameter a. (Right) Dependence of the critical 
dimension A[3, p]/s,p = 1,3 on anisotropy parameter a\ (a 2 = 0) for different values of the 
compressibility parameter a. 

constructed solely of the scalar gradients (see, e.g., Ref. [29] for details) 

F[N,p] = d ll d---d lp e(dM0) n , (52) 

where N = p + 2n is the total number of the fields 9 entering into the operator and p is the 
number of the free vector indices. 

5.2 Composite operators F[N,p]: renormalization and critical di- 
mensions 

Here we shall briefly discuss the renormalization of the composite operators (|52|) which play 
central role in our investigation (complete and detailed discussion of the renormalization of this 
composite operators can be found in Ref. [29J) but first we describe the basis of general theory. 
The necessity of additional renormalization of the composite operators ( 1521) is related to the 
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Figure 6: Dependence of the critical dimension A[3, p]/e,p = 1, 3 on anisotropy parameters «2 
(«! = 0) and «i = a 2 for different values of the compressibility parameter a. 

fact that the coincidence of the field arguments in Green functions containing them leads to ad- 
ditional UV divergences. These divergences must be removed by special kind of renormalization 
procedure which can be found, e.g., in Refs. [201 I2T1 150] . where their renormalization is studied 
in general. The renormalization of composite operators in the models of turbulence is discussed 
in Refs. [23~| [26]. Besides, typically, the composite operators are mixed under renormalization. 

Let F = {F a } be a closed set of composite operators which are mixed only with each 
other in renormalization. Then the renormalization matrix Zp = {Z a p} and the matrix of 
corresponding anomalous dimensions 7^ = {7^/3} for this set are given as follows 

F a = J2 Z ^Ff, lF = Z F l D^Z F . (53) 



Renormalized composite operators are subject to the following RG differential equations 

(V, + (3 g d g + f3 x d x - lv V v )F* = -J2 7apFf, (54) 
which lead to the following matrix of critical dimensions = {A^} 

A F = d k F - A t d F + Y F , A t = -2 + 2e, (55) 
19 




Figure 7: Behavior of the critical dimension A[3, l]/e on anisotropy parameter ot\ and «2 for 
a = (incompressible case) and for a = 0.2. 

where dp a rf'p are diagonal matrices of corresponding canonical dimensions and 7^ is the matrix 
of anomalous dimensions ( 153]) taken at the fixed point. In the end, the critical dimensions of 
the set of operators F = {F a } are given by the eigenvalues of the matrix Ap. The so-called 
"basis" operators that possess definite critical dimensions have the form 

F b a as = E UapFf , (56) 

P 

where the matrix Up — {U a p} is such that A^ = UpApUp 1 is diagonal. 

As was already mentioned, in our case of a scalar admixture with anisotropy and com- 
pressibility the central role is played by the tensor composite operators 8^9 ■ ■ ■ 8i p 9 (8i98i9) n , 
constructed solely of the scalar gradients [29]. It is convenient to deal with the scalar operators 
obtained by contracting the tensors with the appropriate number of the vectors n, namely 

F[N,p] = [(nd)9} p (8 i 98 i 9) n , N = 2n + p. (57) 

Detail analysis shows that the composite operators (15T1) with different iV are not mixed with 
each other in renormalization, that is why the corresponding infinite renormalization matrix 

F[N,p} = Z [NM[N ,, p/] F R [N',p'} (58) 

20 




Figure 8: Dependence of the critical dimension A[4,p]/e,p = 0,2,4 on anisotropy parameters 
«i («2 = 0) and a 2 («i = 0) for different values of the compressibility parameter a (solid lines 
- a = 0, dash lines - a = 0.1, dot lines - a = 0.2, and dash dot lines - a = 0.5). 

is block-triangular, i.e., Z[iv,p][jV'y] = for N' ^ N. Thus, the critical dimensions associated 
with the operator F[N,p] are completely determined by the eigenvalues of the subblocks with 
N' = N. 

In the isotropic well as in the case when large-scale anisotropy is present, the 

elements Z[n, p ] [n, p '] vanish for p < p', therefore the block Z[jy )P ] [jv, p '] is triangular. The same is 
valid for the matrices Up and Ap defined in Eqs. (156 p and fl55l) . Thus, the presence of large- 
scale anisotropy does not affect critical dimensions of the operators (jSTp . Situation is radically 
different in the case when small-scale anisotropy is present, where the operators with different 
values of p mix in renormalization in such a way that the matrix Znv,p] \N,p'] is n °t triangular 
and one can write 

[N/2\ 

F[N,p}= Z [N , pmN -2i]F R [N,N-2l], (59) 

1=0 

where |_7\^/2J means the integer part of the N/2. Therefore, each block of renormalization 
constants with given N is an ( |_7V/2 J + 1) x ( \_N/2\ + 1) matrix. Of course, the matrix of critical 
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Figure 9: (Left) Dependence of the critical dimension A[A,p]/e,p = 0,2,4 on anisotropy pa- 
rameters «i = «2 for different values of the compressibility parameter a. (Right) Dependence 
of the critical dimension A[5,p]/e,p = 1, 3, 5 on anisotropy parameter ot\ (a 2 = 0) for different 
values of the compressibility parameter a. See the caption in Fig. [8] for line identification. 

dimensions fl55l) . whose eigenvalues at IR stable fixed point are the critical dimensions A[N,p] 
of the set of operators F[N,p], has also dimension (\_N/2\ + 1) x ( |_JV/2J +1). 

In what follows, we shall calculate the matrix ^[jv,p][JV',p'l i n one-loop approximation (note 
that in spite of the UV renormalization of our model, where one-loop result is the complete 
solution of the problem because of nonexistence of higher-loop corrections, the renormalization 
of the composite operators has two- and higher- loop corrections). It was analyzed in detail 
in Ref. [22], therefore we shall not repeat it here, and we confine ourself only to the necessary 
information. 

We are interested in iV-th term, denoted as IV [a;; 0], of the expansion of the generating 
functional of one-particle-irreducible Green functions with one composite operator F[N, p] and 
any number of fields 8, denoted as T[x; 8}. It has the following form 

T N (x;8) = —jdxi---j dx N 8(x 1 )---8(x N )(F[N,p](x)8(x 1 )---8(x N )) 1 _ ir (60) 
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Figure 10: Dependence of the critical dimension A[5,p]/e,p = 1, 3, 5 on anisotropy parameters 
a 2 («i = 0) and a\ = a 2 for different values of the compressibility parameter a (see the caption 
in Fig. E]). 

In one-loop approximation the function ( 1601) is given as 

T N = F[N,p} + T^\ (61) 

where is given by the analytical calculation of the diagram in Fig. [3] (see, e.g., Ref. [29] for 
details) and the first term in Eq. (IBTI) represents "tree" approximation. 

After cumbersome but straightforward calculations (see Ref. [29]) one finds the UV divergent 
part of from Eq. (IBTT) in the following appropriate form 

with 

Qi = H Q {A i0 + aA a ) + H x {B iQ + aB a ) , (63) 

where Hj = 2 F\ Q, 1; j + |; — xj ,3 — 0, 1 are corresponding hypergeometric functions, and 
coefficients and Bij for i = 1,2,3,4 and j = 0, 1 are given in Appendix A. 
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Figure 11: Dependence of the critical dimension A[10, p]/e,p = 0,2,4,6 on anisotropy param- 
eter a.\ («2 = 0) and a 2 (a% = 0) for different values of the compressibility parameter a (see 
the caption in Fig. 

The renormalization matrix ^[jv,p] [jvy] is then found from the condition that F^iV, p] from 
Eq. (1611) must be free of UV divergences when are written in renormalized variables, i.e., it does 
not contain poles in e. In the MS scheme one obtains 



gS d A 

6[N,p][N,p+2(i-2)] — <J2i + /r> _ NW ~ 



(64) 



(27r) d 8e 

where i = 1, 2, 3, 4, A = l/(d 2 x 2 (d + 2)(d + 1)), and coefficients Qi,i = 1,2,3,4 are given in 
Eq. (I631) . In the end, from Eq. (164"]) and using definitions fl53|) . one finds the following expressions 
for the elements of the matrix of anomalous dimensions j[N,p] [N',p'] 



7[JV,p][JV,p+2(i-2)] 



(65) 



(2?r) d 4 

for 2 = 1, 2, 3, 4 and the matrix of critical dimensions (155!) has the form 

A[ N ,p]{N,p>] = Ne + 7[iv,p][JV,pT> 

where the asterisk means that 7[jv,p][jv,p'] is taken at the fixed point given by Eqs. fl36l) and ([3j 
Eq. (1661) . which depends on the anisotropy parameters a±,a2, as well as on the compressibility 
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Figure 12: (Left) Dependence of the critical dimension A[10,£>]/£,p = 0,2,4,6 on anisotropy 
parameters «i = «2 for different values of the compressibility parameter a. (Right) Dependence 
of the critical dimension A[ll,p]/e,p = 1, 3, 5, 7 on anisotropy parameter «i («2 = 0) for differ- 
ent values of the compressibility parameter a. See the caption in Fig.[8]for line identification. 



parameter a, is desired one- loop expression for the matrix of critical dimensions of the composite 
operators (IBTj) . 

5.3 Anomalous scaling: one-loop approximation 

The critical dimensions of the operators F[N,p], which we denote as A[N,p], are after all 
equal to the eigenvalues of the matrix of critical dimensions ( 1661) . In the isotropic case or in 
the case with large-scale anisotropy, where the matrix (1661) is triangular, the eigenvalues are 
equal directly to the diagonal elements of the matrix. This fact allows us to assign uniquely 
the concrete critical dimension to the corresponding composite operator even in the case with 
small-scale anisotropy and study their hierarchical structure as a functions of p (see [29J for 
details). In Ref. [29] it was shown that some of the critical dimensions are negative, therefore 
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«2 "l 

Figure 13: Dependence of the critical dimension A[ll, p]/e,p = 1, 3, 5, 7 on anisotropy param- 
eter a 2 («i = 0) and a\ = a 2 for different values of the compressibility parameter a (see the 
caption in Fig. 

they lead to the anomalous scaling (singular behavior of the scaling functions). At the same 
time, the leading role is played by the operators with the most negative critical dimensions 
which are the operators with p = for the structure functions (|44|) with even N, and the 
operators with p = 1 for the structure functions with odd N. 

Thus, the combination of the RG representation with the OPE (1491) leads to the final 
asymptotic expression for the structure functions (jUJ) within the inertial range 

S N (r) = D N/2 r N ^ J2 £{CW (r/L) A ^ + ...}, (67) 

N'<N P 

where p obtains all possible values for given N', Cn>,p are numerical coefficients which are 
functions of the parameters of the model, and dots means contributions by the operators others 
than F[N,p] (see, e.g., [2Tj [29] for details). 

In what follows, we shall investigate the influence of compressibility on the picture found in 
Ref . [29] . Our aim is to find the answer on the question whether the presence of compressibility 
leads to the more pronounced anomalous scaling of structure functions of the scalar field or 
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Figure 14: Behavior of the critical dimension A[10, 0]/e on anisotropy parameters a\ and a-i 
for a = (incompressible case) and for a = 0.2. 

not. Another purpose of the present analysis is to check whether the combine effects defined 
by compressibility and small-scale anisotropy can lead to the more complicated structure of 
critical dimensions than it was studied in Ref. [29] (see also Ref. [?]). One possibility is that the 
pairs of complex conjugate eigenvalues of the matrix of critical dimensions can exist. It leads 
to the oscillation behavior of the corresponding scaling function [29], i.e., the scaling functions 
in Eq. . (1671) would contain terms of the following form 

(r/L) A « {ci cos [A/(r/L)] + c 2 sin [A 7 (r/L)]} , (68) 

where An and Aj are real and imaginary part of A, and c\$ are constants. Another possibility is 
related to the situation if the matrix of critical dimensions cannot be diagonalized and has only 
the Jordan form. Then a logarithmic correction would be involved to the powerlike behavior 
(see Ref. [22]). 

Behavior of the eigenvalues of the matrix of critical dimensions A[N,p] for various values of 
N for different values of the compressibility parameter a and as the functions of the anisotropy 
parameters a.\ and a<i are shown in Figs.HlfT5l It can be immediately seen that only real 
eigenvalues of the corresponding matrix exist as in the incompressible case studied in Ref. [29] . 
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Figure 15: Behavior of the critical dimension A[ll,l]/e on anisotropy parameters a\ and a-i 
for a = (incompressible case) and for a = 0.2. 

At the same time, the hierarchical structure of the critical dimensions for different values of p 
and the same N shown in Ref. [29] is also preserved. The dependence of the critical dimension 
A [2,0] is not shown because it is identically equal to zero. It can be shown either by direct 
calculation or by using the Schwinger equation (see, e.g., Ref. [29]). All the figures show that 
compressibility leads to smaller values of the critical dimensions of the corresponding composite 
operators, thus, at the same time, to the more noticeable anomalous scaling. As it is shown 
in the figures, the effects of compressibility are the most pronounced for small values of the 
anisotropy parameters, as well as for the negative values of these parameters. On the other 
hand, in the limit of large positive values of anisotropy parameters ol\ and the effects of 
compressibility tend to zero. Especially, in the limit ol\^ — > oo the influence of compressibility 
is vanished completely which can be shown by direct analytical investigation. 

6 Conclusion 

In this paper we have analyzed asymptotic behavior of the single-time structure functions S^(r) 
of a passively advected scalar field by the compressible velocity field with small-scale anisotropy 
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by using the field-theoretic renormalization group and the operator product expansion in a 
minimal-subtraction scheme of analytical renormalization. 

It is shown that the leading-order powerlike asymptotic behavior of the single-time structure 
functions of a scalar field within the inertial range in our compressible anisotropic case is given 
by the critical dimensions of the same composite operators as in the incompressible anisotropic 
case but now they acquire rather strong dependence on the compressibility parameter. Further, 
it is shown that when the parameter of compressibility is increasing then the critical dimensions 
of relevant composite operators become smaller, therefore, within our model, we can conclude 
that the anomalous scaling of the structure functions of a passive scalar quantity advected by 
a given velocity field is more pronounced in the compressible stochastic environment than in 
the incompressible one (see Figs.l4l fT5l) . Concrete calculations were done up to the structure 
functions of order N = 11. In our calculations we have not found possible oscillatory modulation 
(related to the possible existence of the complex conjugate eigenvalues of the matrix of critical 
dimensions), as well as we have not found possible logarithmic corrections (related to the fact 
that the corresponding matrix can has only Jordan form) to the leading powerlike asymptotic. 
Thus, all calculated corrections have had purely powerlike behavior. 

It is also shown that the critical dimensions are ordered hierarchically as in the incom- 
pressible case [29], i.e., the compressibility does not disturb the anisotropic structure of the 
critical dimensions rather it shifts them toward smaller values. The largest shift of the critical 
dimensions is present for negative and small positive values of the anisotropy parameters. On 
the other hand, when the anisotropy parameters tend to infinity the effects of compressibility 
vanish. 
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Appendix 

The explicit form of the coefficients A*j and Bij with i = 1, 2, 3, 4 and j = 0,1 (see Eq. (|63|) ) 
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io = p(p-l)d(l + d)((a 1 -a 2 )(2 + d)-(d + 2 + a 2 (l + d) + ai(d 2 -3))x 
+ (d - 2 + a 2 + «i(l + d) + d 2 )x 2 )), 

A n = p (p-l)d(l + d)(2 + d)(l + X )x, 

p{p - l)(-(ai - a 2 )d{l + d)(2 + d) + (1 + d)(2 + d)(a 2 (d - 1) + ttl (d - l) 2 + 
-(d 2 - l)(-2 + cti + a 2 + d + aid + d 2 )x 2 ), 

Bn = p(p-l)(-(l + d)(2 + d)d X -(d + 2)(d 2 -l)x 2 ), 

. 20 = rf(4n(n-l)(l + x)(3(2 + rf)(x-a 1 ) + a 1 (rf 2 -l)x 

+3a 2 (2 + d + x (l + d))) - (1 + <0(p(p - !)(- rf ( 2 + + «i(d(2 + d) 

+2 X (1 + x)) - « 2 (1 + X)K2 + d(l + X )) - 2 X )) 

-2((2 + d)x(a 2 - ai + (a 2 + a ± (d - 1) + d)x) 

-2(-(2 + d)x((d - 3)x - 3) + a 2 (l + x)(6 + 3d + x + 2d X ) 

+ax(d 2 x + (X - 6)(1 + x) - d(3 + x(2 + x))))p)n)), 

A 21 = d(-12n(n-l)(l + x)(2 + d)x + (l + d)(-p(p-l)d(2 + d)x 
+ 12(l + x)(2 + d)xpn)), 

30 = _4 n ( n _i)(_3( ai _ a2 )rf(2 + d) + (2 + d)(3d-3« 2 + 6a 2 d 
+ai(3 + (d - 5)d))x + (7d - 6 - 3a 2 + (3 + 3a 2 - d)d 2 
+a 1 (d - 1) 2 (1 + d))x 2 ) + (1 + d)(p(p - l)((ai - a 2 )d 2 (2 + d) 
-d(2 + d)(2« 2 (-l + d) + d)x - (d - l)(-2(ai + a 2 ) + « 2 d 2 )x 2 ) 
-2((2 + d)x((a 2 - «i)d + (a 2 (d - 1) + a ± (d - l) 2 + d)x) 
+2(3(«i - a 2 )d(2 + d) - (2 + d)(a 1 (d - 3)(d - 1) + 3d - 3a 2 + 5a 2 d)x 
+ (d - l)(-6 -ai-a 2 + (an - 2a 2 - l)d + d 2 )x 2 )p)n), 
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-12ra(n - l)(-(2 + d)d X + (2 - 3d - 2d 2 ) X 2 ) - (1 + d)(-p(p - l)d 2 (2 + d) X 
+2(-(2 + d)d X 2 + 6(2 + d) X (d + (d - l)x)p)n), 

= -2d(2 + d)(2n(n - l)(-6«i(2 + d) + «i((d - 5)d - 12)x 

+6a 2 (l + X )(2 + d + x + d X ) + *(12(1 + x) + <*(6 + 5 X - d X ))) 

+ (1 + d)(dx(-ai + a 2 + x + « 2 x) - 2(-a 1 (2 + d + 2 X ) + x(2 + d + 2*) 

+a 2 (l + x)(2 + d + dx))p)n), 

2d(2 + d)(12n(n - 1)(1 + x)(2 + d) X - (1 + d)(-d X 2 + 2( X (2 + d + 2 X ))p)n), 

= 2(2 + d)(2n(n - l)(-6(«i - a 2 )d{2 + d) + (6(2 + d)(d + a 2 (2d - 1)) 
+«i(12 + d(-12 + (d - 5)d)))x + (-12 + d(12 - (d - 5)d) 
+a 2 (d - 6 + 7d62))x 2 ) + (1 + d)(dx(-«id + d X + a 2 (d + 2(d - l)x)) 
-2(-(«i - a 2 )d(2 + d) + (-2«i(d - 1) + d(2 + d) + 2a 2 (d - 1 + d 2 )) X 
+ (d-l)(2 + a 2 d) X 2 )p)n), 

B 31 = 2(2 + d)(12n(n-l)(-x(d(2 + d)(l + x)-2x)) 
+ (1 + d)(-dV + 2(xK2 + d) + 2(d - l)x))p)n), 

40 = 4n(n-l)d(2 + d)(a 2 (4 + d)(l + x)(2 + d + x + rfx) 

-ai(8(l + x) + d(G + d + 5x)) + x(8(l + x) + d(G + d + 5 X ))), 

An = -4n(n-l)d(2 + d)x(8(l + x) + c?(6 + d + 5x)), 

An(n - 1)(2 + d)((«i - a 2 )d(2 + d)(4 + d) - (2 + d)(e*i(4 - 5d) 
+ (4 + d)(d + a 2 (2d - l)))x - ((2 + d)(5d - 4) + a 2 (-4 + d 2 (5 + d)))x 2 ), 

B 41 = 4n(n- l)(2 + d)(d(2 + d)(4 + d)x+ (-8 + 6d + 5d 2 )x 2 ), 
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